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The precise physical mechanism describing the emergence of the seeds of cosmic structure from a 
perfect isotropic and homogeneous universe has not been fully explained by the standard version of 
inflationary models. To handle this shortcoming, D. Sudarsky and collaborators have developed a 
proposal: the self-induced collapse hypothesis. In this scheme, the objective collapse of the inflaton 
wave function is responsible for the emergence of inhomogeneity and anisotropy at all scales. In 
previous papers, the proposal was developed with an almost exact de Sitter space-time approxima¬ 
tion for the background that led to a perfect scale-invariant power spectrum. In the present article, 
we consider a full quasi-de Sitter expansion and calculate the primordial power spectrum for three 
different choices of the self-induced collapse. The consideration of a quasi-de Sitter background allow 
us to distinguish departures from an exact scale-invariant power spectrum that are due to the inclu¬ 
sion of the collapse hypothesis. These deviations are also different from the prediction of standard 
inflationary models with running spectral index. Comparison with the primordial power spectrum 
and the CMB temperature fluctuation spectrum preferred by the latest observational data is also 
discussed. From the analysis performed in this work, it follows that most of the collapse schemes 
analyzed in this paper are viable candidates to explain present observations of the CMB fluctuation 
spectrum. 


I. INTRODUCTION 

Recent observations of the Cosmic Microwave Back¬ 
ground (CMB) radiation are one of the most powerful 
tools to study the early universe and also to provide a 
method to settle the value of the cosmological parame¬ 
ters. In the last 20 years, there has been a lot of improve¬ 
ments in the measurements precision of the CMB radi¬ 
ation anisotropies. Furthermore, the agreement between 
theory and observations has strengthened the theoretical 
status of inflationary scenarios among cosmologists. 

In the standard inflationary paradigm, the emergence 
of all structures in our universe like galaxies and galaxy 
clusters are described by a featureless stage represented 
by a background Friedmann-Robertson-Walker (FRW) 
cosmology with a nearly exponential expansion driven by 
the potential of a single scalar field and from its quan¬ 
tum fluctuations characterized by a simple vacuum state. 
However, when this scenario is considered more carefully, 
a conceptual problem emerges regarding a change in the 
initial symmetries of the universe. This is, from a highly 
homogeneous and isotropic initial state that characterizes 
the quantum perturbations of both the classical back¬ 
ground inflaton and space-time, the universe ends in a 
state with “real” inhomogeneities and anisotropies. In 
other words, if one considers quantum mechanics as a 
fundamental theory, then it is appropriate to use it in 
order to describe the universe as a whole; therefore, ev¬ 
ery classical description of the universe shall be associ¬ 
ated to an imprecise highly complex quantum mechan¬ 


ical state. Moreover, the observed universe, at certain 
scales, does exhibit inhomogeneous and anisotropic fea¬ 
tures; consequently, its quantum description in terms of 
a quantum state must encode these non-symmetrical as¬ 
pects. On the other hand, the dynamics of the standard 
inflationary paradigm does not contain any aspect that 
can be responsible for breaking the initial symmetries of 
the early quantum state, which happened to be perfectly 
homogeneous and isotropic. In this sense we consider 
that the standard inflationary paradigm is incomplete. 
D. Sudarsky and collaborators (ll-flij have developed one 
proposal to handle these shortcomings. The way to deal 
with the problem is introducing a new ingredient into the 
inflationary account of the origin of cosmic seeds: the 
self-induced collapse hypothesis , i.e. a scheme in which 
an internally induced collapse of the wave function of the 
inflaton field is the mechanism by which inhomogeneities 
and anisotropies arise at each particular length scale. 

The collapse proposal was inspired by previous ideas 
of R. Penrose and L. Diosi [llj 14], which regarded the 
collapse of the wave function as an actual physical pro¬ 
cess induced by gravity (instead of just an artifact of the 
description of Physics). At this point of the discussion, 
we do not know exactly what kind of physical mechanism 
would lie behind what, at the semiclassical level we are 
working, looks like a spontaneous collapse of the wave 
function. We assume that the effect is caused by un¬ 
known quantum aspects of gravitation. Essentially, the 
collapse hypothesis simply sustains that something in¬ 
trinsic to the system, i.e. independent of external agents 
(e.g. observers), induces the collapse or reduction of the 
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quantum mechanical state of the system. Various propos¬ 
als of that sort have been developed EMI, and might 
well be compatible with the self-induced collapse of the 
inflaton’s wave function that we are considering. The 
proposal is, at this stage of the analysis, a purely phe¬ 
nomenological scheme. It does not attempt to explain 
the process in terms of some specific new physical the¬ 
ory, but merely gives a rather general parametrization of 
the quantum transition involved. We will refer to this 
phenomenological model as the collapse scheme. 

Here, it is worthwhile to mention that the previous 
conceptual problem is sometimes known in the literature 
as the quantum-to-classical transition of the primordial 
perturbations. As a matter of fact, a partially under¬ 
standing of such issue has been gained by relying on the 
decoherence framework and the fact that the initial vac¬ 
uum state of the inflaton evolves into a highly squeezed 
state EMI; in particular, it is usually argued that the 
predictions from the quantum theory, characterizing the 
inflaton fluctuations, are indistinguishable from those of 
a theory in which the random fluctuations are the result 
of a classical stochastic process [20|. However, this argu¬ 
ment alone cannot explain the fact that a single (classical 
and random) outcome emerges from the quantum theory 
[2lJ. In other words, decoherence (and the squeezing of a 
quantum state) cannot solve the quantum measurement 
problem [22l 1231 ], which in the cosmological setting is am¬ 
plified, i.e. there is no clear way to define entities such as 
observers, measurement devices, environmental degrees 
of freedom, etc. Other cosmologists seem to adopt the 
Everett “many-worlds” interpretation of quantum me¬ 
chanics when confronted with the quantum-to-classical 
transition in the inflationary universe (although a de¬ 
tailed and precise formulation of such posture is still not 
well represented in the literature). In the Everettian for¬ 
mulation, reality is made of a connected weave of ever- 
splitting worlds, each one realizing one of the alterna¬ 
tives that is opened by what we would call a quantum- 
mechanical measurement. Regarding this point, we want 
to mention that Everett’s interpretation has not com¬ 
pletely solved the measurement problem yet. In fact, 
there is a mapping between what in that approach would 
be called the splittings of the worlds, and what would be 
called the “measurements” in the Copenhagen interpre¬ 
tation. Therefore, every question that can be made in 
the orthodox interpretation has a corresponding one in 
the Everettian one. That is, the specific issues regarding 
the measurement problem would be: When and why does 
a world splitting occur?, under what circumstances does 
it occur? What constitutes a trigger? Furthermore, even 
if one could bypass those conundrums, a precise deriva¬ 
tion of the Born rule (a crucial aspect to connect the 
theory’s predictions with the experiments) and a clear 
justification for the choice of a particular basis in which 
the splitting takes place is unknown within Everett’s in¬ 
terpretation (24l [25|. 

Given the previous discussion, an objective reduction 
of the quantum state, characterizing the inflaton field, 


seems to be a plausible option for addressing the problem 
at hand. The detailed analysis of the original proposal 
at the conceptual and technical level can be consulted in 
Refs. 00. 

Furthermore, even though there are well-known models 
characterizing the self-induced collapse of the wave func¬ 
tion E103 in a generic (non-cosmological) context, the 
relativistic framework for the objective collapse models 
is still under development [26|, [23 ■ In this work, we will 
follow a more pragmatical approach and characterize the 
post-collapse quantum state by the expectation values 
of the field and its momentum conjugate, without rely¬ 
ing on some particular collapse mechanism. On the other 
hand, there are still various possibilities regarding the de¬ 
scription of the quantum expectation values in the post¬ 
collapse state, we will refer to these various prescriptions 
as collapse schemes. In a first attempt E, 3j, two differ¬ 
ent schemes were considered: one in which, after the col¬ 
lapse, both expectation values of the field and momentum 
variable are randomly distributed within their respective 
ranges of uncertainties in the pre-collapsed state, and 
another one in which it is only the conjugate momentum 
that changes its expectation value from zero to a value 
in its corresponding range, as a result of the collapse. In 
later works 0, f|, another scheme was considered moti¬ 
vated by the correlation between the field variable and its 
conjugated momentum existing in the pre-collapse state; 
in this scheme the collapse is characterized in terms of 
the Wigner function. 

In previous papers 00 the collapse proposal was de¬ 
veloped using an “almost” exact de Sitter expansion fac¬ 
tor during the inflationary period [i.e., it was assumed 
that the (cosmic) time derivative of the Hubble factor H 
was exactly a constant]. In consequence, the primordial 
power spectrum, modified by the collapse hypothesis, re¬ 
sulted in a form P(k ) = A s C{k) 1 this is, a power spec¬ 
trum with scalar spectral index n s = 1, and with C(k ) 
a function representing the modification induced by the 
collapse hypothesis; in particular, the function C(k) de¬ 
pends on the time of collapse of each mode of the infla¬ 
ton field. Furthermore, the temperature and polarization 
power spectrum of the CMB were also modified, there¬ 
fore, the proposal predictions could be tested against cur¬ 
rent observational data 0 . In such way, it was possible 
to recover an exact scale-invariant power spectrum for 
some values of the collapse parameters that are related 
to the time of collapse of each mode of the inflaton field. 

However, recent data reported by Planck [20 and 
WMAP [20], constrain the value of the spectral scalar 
index to n s = 0.9603 ± 0.0073. Therefore, in this pa¬ 
per, we go one step beyond and calculate the primor¬ 
dial power spectrum for different collapse schemes in a 
full quasi-de Sitter background, i.e. by considering small 
time variations of H. In this manner, we will obtain an 
expression of the form P(k) = A s k n ’ > ~ 1 Q(k)- 1 naturally, 
with n s 1, and also with Q(k) a function introduced 
by the collapse hypothesis, which can be reduced to the 
conventional phenomenological expression for some val- 
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ues of the collapse parameters. We emphasize that this 
work is not only a matter of technical improvement, but 
also helps to separate the features in the power spec¬ 
trum that can be attributed to the collapse of the wave 
function and the aspects that are only due to consider¬ 
ing a quasi-de Sitter background. That is, in previous 
works the prediction of the form P(k) = A s C(k) does 
not allow one to exactly identify the dependence on k at¬ 
tributed to the collapse hypothesis and to the spectral in¬ 
dex; meanwhile, in the prediction obtained in the present 
manuscript, P(k) = A s k ne ~ 1 Q(k), one can plainly rec¬ 
ognize and separate the two kind of dependencies. We 
think that this last feature is important since, at the phe¬ 
nomenological level, one wishes to learn more about the 
unknown collapse mechanism, so it is of great significance 
to clearly identify its characteristics within the theoreti¬ 
cal predictions that will be confronted with observational 
data. 

Additionally, our prediction for the power spectrum 
allows departures from the traditional inflationary ap¬ 
proach that can be tested observationally. Moreover, as 
we will show later in the paper, since our model is con¬ 
ceptually different from the standard one, we cannot fol¬ 
low the traditional method of simply evaluate the power 
spectrum, obtained from pure de Sitter inflation, at the 
so-called “horizon crossing” and, in this way, achieve a 
power spectrum that is equivalent in shape to the cor¬ 
responding one calculated in a quasi-de Sitter stage. In 
order for our model to yield a conceptually consistent pre¬ 
diction, we must perform the full calculation in a quasi-de 
Sitter background. 

On the other hand, in a recent work [9f we have an¬ 
alyzed the collapse hypothesis in the case that the col¬ 
lapse occurs during the radiation epoch; and we have 
shown that if one considers an almost exact de Sitter ex¬ 
pansion for the inflationary era, then the model cannot 
account for the present observational CMB temperature 
fluctuation spectrum. This last statement also applies 
to the case where the collapse happens during inflation 
(which is the case of the present paper). In fact, a sim¬ 
ple calculation of the y 2 value for the models presented 
in Ref. [7j, where an almost exact de Sitter background 
was considered using the CMB temperature data released 
by Planck 28(, the CMB polarization data reported by 
WMAP [2S |, CMB temperature data for high values of l 
reported by the Atacama Cosmology Project (ACT) {30] 
and the South Pole Telescope (SPT) [3l|, yields a value 
that is several orders of magnitude above the expected 
reasonable value. The difference between the analysis 
performed in the former paper and this simple calcula¬ 
tion is that the data obtained by the Planck, ACT and 
SPT collaborations are much more accurate for small an¬ 
gular scales than the 7-years release of WMAP. 

The paper is organized as follows: In Sec. UD we re¬ 
view the senriclassical gravity approach, in which only 
the perturbations of the inflaton field are quantized and 
obtain the corresponding linearized Einstein’s equations. 
In Sec. mu we perform the quantization of the infla¬ 


ton field in a quasi-de Sitter background. Furthermore, 
we relate the initial curvature perturbation with the pa¬ 
rameters characterizing each collapse scheme; we con¬ 
sider three different choices for the quantum collapse. In 
Sec. EO we relate the CMB observational quantities with 
the primordial spectrum modified with the collapse hy¬ 
pothesis in a quasi-de Sitter background. In Sec. El we 
compare the primordial power spectrum obtained in Sec. 
II VI with the phenomenological expression in standard in¬ 
flationary models. In Sec. ED we plot the primordial 
power spectrum obtained in this paper for some partic¬ 
ular values of the collapse parameters and compare it 
with the primordial spectrum preferred by the data. In 
Sect. IVIII we present the prediction for the CMB tem¬ 
perature fluctuation spectrum and show that there are 
different predictions corresponding to the three collapse 
schemes proposed along with distinct values for the time 
of collapse. Finally, in Sec. IVIIII we summarize the main 
results of the paper and present the conclusions. 

Regarding notation and conventions, we will work with 
signature (—,+,+,+) for the metric; primes over func¬ 
tions will denote derivatives with respect to the confor¬ 
mal time 77 , and we will use units where c = h = 1 but 
keep the gravitational constant G. 


II. SEMICLASSICAL GRAVITY AND 
LINEARIZED EINSTEIN’S EQUATIONS 

The purpose of this section is to present our view re¬ 
garding the relation between the space-time description 
in terms of the metric and the quantum degrees of free¬ 
dom (DOF) of the matter fields, represented by the in¬ 
flaton. First, we will introduce the physical point of view 
of such relation and then derive the corresponding equa¬ 
tions. 

We have proposed that the mechanism by which pri¬ 
mordial anisotropies and inhomogeneities arise is a self- 
induced collapse of the inflaton wave function. As a con¬ 
sequence, the post-collapse quantum state must not be 
an homogeneous and isotropic state, this is, it is not an 
eigen-state of the linear and angular momentum oper¬ 
ators. One could then assume that the post-collapse 
state was obtained from a particular collapse mechanism, 
and then compute the corresponding observables in that 
state. The question now would be: In the context of the 
inflationary scenario, what are the appropriate observ¬ 
ables that result from the quantum theory? 

One possible approach would be to assume that both 
metric and matter perturbations are described by a 
quantum field theory constructed on a classical unper¬ 
turbed background; in the inflationary universe, this ap¬ 
proach corresponds to the quantization of the so-called 
Mukhanov-Sasaki variable. Henceforth, if one assumes 
a particular collapse mechanism, which somehow modi¬ 
fies the standard unitary evolution given by Schrodinger 
equation, then the dynamics of the Mukhanov-Sasaki 
variable would induce non-standard predictions for the 
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observational quantities (e.g. the spectrum of the tem¬ 
perature anisotropies). This scheme was developed in 
Refs. [3214351 ] for the inflationary universe. 

Another potential prescription to relate the quantum 
DOF with the observational quantities, is to rely on the 
semiclassical gravity approximation represented by Ein¬ 
stein semiclassical equations G a b = 87 rG(T a b); within this 
framework, the metric is described in a classical way, 
while the matter DOF are modeled by a quantum field 
theory in a curved classical background. Then, during 
inflation, the semiclassical equations enable to relate the 
inflaton quantum perturbations with the corresponding 
ones from the classical metric. However, assuming a 
particular collapse mechanism, which again can be en¬ 
visioned as a modification of the Schrodinger equation, 
would not alter the evolution of the metric perturbation; 
indeed, the dynamics of the modes 

characterizing the quantum field, representing the in¬ 
flaton, would be modified, but the metric perturbation is 
always a classical object, and, thus, its evolution is not 
dictated by the modified Schrodinger equation. Assum¬ 
ing a particular collapse mechanism, would only modify 
the initial conditions of the motion equation for the met¬ 
ric perturbation, which again is always described at the 
classical level; in the context of inflation, this was ana¬ 
lyzed in Ref. 0 (by assuming pure de Sitter inflation). 
Nevertheless, the initial condition for the motion equa¬ 
tion of the metric perturbation, will contain the informa¬ 
tion regarding that a collapse of the wave function of the 
inflaton has occurred. 

One main advantage of using the semiclassical picture 
is that the description and treatment of the metric (both 
the background and its perturbations) remains “classi¬ 
cal” at all times. As a consequence, there is no issue 
with the “quantum-to-classical transition” in the sense 
that one needs to justify going from “metric operators” 
(e.g. T) to classical metric variables (such as f). The 
fact that the space-time remains classical is particularly 
important in the context of models involving dynamical 
reduction of the wave function, as such collapse or reduc¬ 
tion is regarded as a physical process taking place in time 
and, therefore, it is clear that a setting allowing consid¬ 
eration of full space-time notions is preferred over, say, 
the timeless settings usually encountered in canonical ap¬ 
proaches to quantum gravity (see Ref. [36] for a complete 
review on the problem of time in quantum gravity). 

Another aspect of the semiclassical approximation, is 
that it is not valid during the collapse. The reason is that, 
as is well-known, introducing a dynamical collapse gener- 
ically violates the conservation of energy, so the diver¬ 
gence of the energy-momentum tensor does not vanish, 
V a (T' ab ) 7 ^ 0. If the divergence of the energy-momentum 
tensor does not vanish, and it is equated to the Einstein 
tensor, then of course the latter’s divergence does not 
vanish either, V a G ab 7 ^ 0, which evidently is a prob¬ 
lem since we know that the divergence of the Einstein 
tensor must be zero. Therefore, during the collapse, we 
cannot say how the modified dynamics of the quantum 


fields, provided by a collapse mechanism, affects the clas¬ 
sical metric perturbations that are directly related to the 
observables. However, the semiclassical gravity approxi¬ 
mation is valid before and after the collapse, which cor¬ 
respond to the cases of interest for the present work. 

More specifically, before the collapse, the initial state 
of the universe, characterized by a few e-folds after in¬ 
flation has begun, is described by both the homogeneous 
and isotropic (H&I) Bunch-Davies vacuum and the H&I 
Friedmann-Robertson-Walker space-time. Afterwards, 
at some point during the inflationary epoch, the quantum 
state of the matter fields spontaneously changes to a new 
quantum state (i.e. the post-collapse state); this change 
is induced by some physical (but unknown) mechanism. 
However, the state resulting from the collapse needs not 
to share the same symmetries as the initial state. After 
the collapse, the gravitational DOF are assumed to be, 
once more, accurately described by Einstein’s semiclas¬ 
sical equations. Nevertheless, ( T a b ) evaluated in the new 
state does not generically posses the symmetries of the 
pre-collapse state; hence, we are led to a new geometry 
that is no longer H&I. We have presented just a very 
brief summary of the semiclassical picture and its rela¬ 
tion with the collapse hypothesis during inflation. The 
full formalism has been developed in Ref. j5] and we in¬ 
vite the reader to consult such reference. We also should 
mention that we are not advocating that semiclassical 
gravity must be regarded as a fundamental theory, we 
are using it as an appropriate approximation given the 
energy scales of the inflationary universe (~ 10 16 GeV), 
i.e. we are treating it as a quantum field theory in a 
curved space-time. 

It is also worthwhile to mention a few remarks regard¬ 
ing the tensor modes and the semiclassical gravity ap¬ 
proach. Last year, the BICEP2 collaboration reported a 
detection of the primordial B-mode polarization consis¬ 
tent with the prediction of standard inflationary models 
i~37i . However, other authors pointed out that without 
an accurate dust map, it is not possible to discern be¬ 
tween dust polarization and polarization due to primor¬ 
dial gravity waves |38l440[ . More recently, the Planck col¬ 
laboration analysis of the angular power spectrum of po¬ 
larized dust emission at intermediate and high Galactic 
latitudes [4l| showed that the dust extrapolated power 
spectrum (obtained by extrapolating Planck 353 GHz 
data to 150 GHz) is of the same magnitude as the B- 
mode polarization power spectrum reported by the BI¬ 
CEP collaboration. Finally, a recent joint analysis of the 
BICEP2/Keck and Planck collaborations shows no evi¬ 
dence for primordial B-mode polarization at low /, mean¬ 
while for high /, they have found evidence for B-modes 
that are originated by gravitational lensing [42j . 

On the other hand, in our approach, the source of 
the curvature perturbations lies in the quantum inhomo¬ 
geneities of the inflaton field (after the collapse). Once 
the collapse has taken place, the inhomogeneities of the 
inflaton feed into the gravitational DOF leading to per¬ 
turbations in the metric components, in particular is a 
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direct source of the scalar perturbations. However, the 
metric itself is not a source of the self-induced collapse. 
Therefore, as the scalar field does not act as a source 
for the metric tensor modes, at least not at first order 
considered here, the analysis concerning the amplitude 
of the primordial gravitational waves should be done at 
the second order in the perturbations; such analysis is be¬ 
yond the scope of this paper and would be the subject of 
future research. Furthermore, in Ref. |43| . the semiclas- 
sical gravity approximation plus a collapse of the infla- 
ton’s wave function results in an undetectable amplitude 
for the primordial gravitational waves; however, those 
authors consider that the state collapses on a spacelike 
hypersurface for all wavelengths modes, this contrasts 
with our view in which the time of collapse depends on 
the mode’s wavelength. 

Now that we have established the conceptual relation 
between the matter and geometry fields, we will proceed 
to find the main equation which will illustrate this con¬ 
nection under the collapse hypothesis. 

In the inflationary regime, the dominant type of matter 
is modeled by a scalar field <j> called the inflaton with 
a potential V responsible for the accelerating expansion 
of the universe. At the end of the inflationary epoch, 
the universe follows the standard Big Bang evolution, 
the transition mechanism is provided by the reheating 
period. 

The inflationary universe is modeled by the action of 
a scalar field minimally coupled to gravity: 


sift, gab] = J d A Xy/~^g 


I 67 tG 


R [g] 


- -V a (j)\/ b (j)g ab - V[4>\ 


(1) 


Varying Eq. © with respect to the metric yields the 
field equations G a b — 8nGT a b. 

We proceed the analysis, in the standard fashion, sep¬ 
arating the metric and the scalar field into a background 
(which is perfectly homogeneous and isotropic) plus a 
small perturbation, i.e. g ab = g^b + dg a b and, using con¬ 
formal coordinates, <f> = 4>o(v) + <5</>(x, rf). 

Focusing first on the background, the space-time is 
characterized by a flat FRW space-time. Einstein equa¬ 
tions for the background Gqq = 87 tGTq^ = 8irGa 2 p 
yield Friedmann equations: 3"H 2 = 87 rGa 2 p, where B = 
a'{rf) / a{rf). The major contribution to the energy density 
p comes from the inflaton potential V. 

In the slow-roll inflationary model, the conformal Hub¬ 
ble factor is characterized by B ~ — 1 / [77 (1 — eh)], with 
eh = 1 — B'/B 2 the Hubble slow-roll parameter, which 
during inflation 1 eh — constant. Note that in Ref. 
QJ, B = — 1/??, this is, the background space-time is 
strictly de Sitter and leads to a final spectrum with 
n s = 1 ; this is different from the present paper where 
we will be considering a quasi-de Sitter background and 


that will lead us to a value for the scalar spectral index 
n s 7 ^ 1 . 

Furthermore, since we will work with the slow-roll 
approximation, then the motion equation for the back¬ 
ground field can be approximated by 2>B4>'q = —a^d^V. 
Additionally, it is convenient to introduce the potential 
slow-roll (PSR) parameters 


_ Ml (d^V\ 
^ ~ 2 \V ) 


Sy = Mp 



( 2 ) 


Thus, by assuming Ey,5y -C 1, one identifies the re¬ 
gion in the potential where the slow-roll approximation 
is valid. Furthermore, during slow-roll inflation eh — ey. 

The normalization of the scale factor will be set to 
a = 1 at the present cosmological time. The inflationary 
era would come to an end at 77 = r] r « —10 -22 Mpc, that 
is, the conformal time during the inflationary era is in 
the range —oo<rj<ri r . 

Next we focus on the perturbations. It is usual to 
decompose the metric fluctuations in terms of its scalar, 
vector and tensor components. In the case of interest 
for this article, we concern ourselves only with the scalar 
perturbations. The scalar metric perturbations in a FRW 
background space-time are generically described by the 
line element: 

ds 2 = a 2 (r)){ — (1 — 2 ip)drj 2 + 2(diB)dx l dr] + 

+ [(1 — 2 ip)Sij + 2didjE]dx l dx j }. (3) 

Since we will work in the semiclassical framework it 
is convenient to work with the gauge-invariant quanti¬ 
ties known as the Bardeen potentials, defined as $ = 
ip + ±[a(B - E')}' and <F = ^ + B{E' - B). In a sim¬ 
ilar way, the perturbations of the inflaton can be mod¬ 
eled by the gauge-invariant fluctuation of the scalar field 
( 77 , x) = S(j> + (p ' 0 {B - E'). 

In I© is shown that combining the perturbed Einstein 
equations (in the absence of anisotropic stress) and the 
slow-roll motion equation, one obtains: 


V 2 ^ + p^ = 4ttG^^ ,( gi) , (4) 

where p = Ti 2 — . In Fourier space, Eq. © reads 


= VT u p{ i - „) aSUri)lGr ’• (5) 

with H the Hubble factor and Ml = 1 / 87 rG the reduced 
Planck’s mass; also, we have used the definition of Ey. 
Friedmann equation and the slow-roll approximation for 
<Pq. Additionally, from the definition of p, one can check 
that p = ehE. 2 . During most of the inflationary regime, 
the inequality k 2 p is satisfied (both when \kg\ 1 
and | A; 7 /1 <C 1), it gets violated when eh starts depart¬ 
ing from being a constant an approaching unity; in other 
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words when inflation is coming to an end. However, since 
modes of observational interest are bigger than the Hub¬ 
ble radius (\krj\ <C 1 ) while the inflationary phase is still 
going on, the approximation k 2 /r remains valid. Thus, 
Eq, (0 is approximated by 


'I'kO?) - kO?) (GI) - ( 6 ) 

Finally, since we mentioned that we will rely on the semi- 
classical gravity framework, Eq. 0 can be generalized 
to 


’M??) ^ \[^fj^~k2 a (S^' k {v) (G1) )- (?) 

Equation © is the main result of this section; the dif¬ 
ference with respect to the perfect de Sitter case will be 
reflected in the motion equation for 6(j)' k (r])( GI \ and also, 
in the fact that H and €h are strictly not constant. 

Let us note that Eq. © is expressed in terms of gauge- 
invariant quantities 'I'k and 5(f>' k (rj)( GI \ Nevertheless, in 
the longitudinal gauge, T represents the curvature per¬ 
turbation, and is related to <5</> in the exact same way 
as in Eq. © 0 . Thus, expression © serves to illus¬ 
trate what we mentioned at the beginning of the section, 
namely, that the quantum treatment is all in the matter 
fields, which during inflation is dominated by the infla- 
ton, while the curvature perturbation is always a classical 
quantity. 


we can expand the action (1TJ up to second order in the 
rescaled variable (i.e. up to second order in the scalar 
field fluctuations) 



The canonical momentum conjugated to y is n = 
dSC^/dy' = y' — ( a'/a)y = a5(j)'. The field and mo¬ 
mentum variables are promoted to operators satisfying 
the equal time commutator relations [y(x,77), 7 r(x', 77)] = 
hl(x-x') and [y(x, 77), y(x', 77)] = [tt(x, rj), 7 r(x', 77)] = 0. 
Expanding the field operator in Fourier modes yields 


y(v, x ) = jjtYl ^( ? ?) eIkx > ( 9 ) 

k 

and analogous expression for 71 ( 77 , x), note that the sum 
is over the wave vectors k satisfying kiL = 2irrii for i = 
1,2,3 with rii integer and y k (r]) = y k {i])a k + yZ(r})a'_ k 

and 7 r k ( 77 ) = 3 fc(??)ak + 9 k (y)ai k , with g k (r/) = y' k {rj) - 
TLyk{v)- From the previous expression it is clear that 
we are taking the quantization on a finite cubic box of 
length L, at the end of the calculations we will go to the 
continuum limit ( L —> 00 , k —> cont.). The equation of 
motion for yk{if) derived from action 0 is 


III. QUANTUM THEORY OF FLUCTUATIONS, 
COLLAPSE SCHEMES AND THE PRIMORDIAL 
CURVATURE PERTURBATION 

In this section, we will present the quantum theory of 
the field variables and characterize the collapse proposal. 
As mentioned in the introduction, we will focus on the 
pragmatical approach first exposed in lj; more on this 
pragmatical approach can be consulted in Refs. [3, 4HU ■ 
Next, we will introduce the collapse schemes and finally 
we will find the expression for the curvature perturbation 
for the three schemes considered. 


y'k(ii) + ( fc2 _ ~ + a 2 d l* v ) = °- ( 10 ) 

Since H = — 1 / [77 (1 — £//)], one finds that a"/a ~ 
(2 + 3 e#)/? 7 2 ; additionally, using the definition of Sy, 
Friedmann’s equation and the explicit form of 7i , one 
has a 2 <9^U ~ 36v/v 2 - Therefore, the motion equation 
Eq. (flUl) is rewritten as 

y'kiv) + (k 2 - 2 + — —^yk(v) = 0 - ( 11 ) 


A. Quantum theory of perturbations 

Given that we are working within the semiclassical 
gravity framework, in which only the matter fields are 
quantized, and that the self-induced collapse generates 
the curvature perturbation, our fundamental quantum 
variable will be the fluctuation of the inflaton field 
<5</>(x, 77 ); thus, we will consider the quantum theory of 
<5<(>(x, 77 ) in a curved background described by a quasi-de 
Sitter space-time. Furthermore, it will be easier to work 
with the rescaled field variable y = a5(p ; consequently, 


The solution of Eq. m is fixed by the canoni¬ 
cal commutation relations between y and if, which give 
[&k, = T 3 6 k .k', thus y k {rj) must satisfy y k g* k -ylg k = 

i for all k at some time 77 . The choice of y k (y) corresponds 
to the choice of a vacuum state for the field, which in the 
present case would be the so-called Buncli-Davies vac¬ 
uum characterized by 


y k (y) 


-7TT] 


1/2 

e »[^+l/ 2 ] Or/ 2 ) ^ ( 1 ) ( _ fe?7) , 


( 12 ) 


4 
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where v = 3/2 + eh ~ Sv and H^\—krj) is the Han- 
kel function of the first kind of order z/Q The solution 
involves a phase e *[''+ 1 / 2 K 7r / 2 ) that we will drop from cal¬ 
culations as it has no observational consequence. 

We note that in the case of an exact de Sitter universe, 
the motion equation would correspond to setting eh = 
5v = 0 in Eq. ([Till . and indeed that was the case studied 
in Refs. |lj, y|. The fact that the motion equation HID 
now involves the slow-roll parameters will lead, as we 
will show in the rest of this article, to a prediction for 
the scalar spectral index that is generically n s ^ 1. 

The self-induced collapse hypothesis is based on con¬ 
sidering that the collapse acts similar to a “measure¬ 
ment” (clearly, there is no external observer or detec¬ 
tor involved). This lead us to consider Hermitian opera¬ 
tors, which in ordinary quantum mechanics are the ones 
susceptible of direct measurement. Therefore, we sepa¬ 
rate y^{rf) and dk (rj) into their real and imaginary parts 
= yk R {y)+iy^ 1 (v) and 7 T k (? 7 ) = ^ R {ri)+iW(jl) in 
this way the operators y^' 1 (y) and d k R (y) are Hermitian 
operators. Thus, 


= V / 2Re[ 2 / fc (7 ? )a k / ], 

(13a) 

n RI (r)) = V2Re{g k (r))a RJ ], 

(13b) 


where = (fik + a_k)/-\/ 2 , a£ = —i(dk — d_k)/-\/ 2 - 

The commutation relations for the a k ’ are non¬ 
standard 

[o£' J ,a£; Jt ] = l 3 * * * (4,k' ± 4,-k'), (14) 

where the + and the — sign corresponds to the com¬ 
mutator with the R and I labels respectively; all other 
commutators vanish. 

Up to this point, we have proceeded in a similar way 
to the traditional one, except that we are treating at the 
quantum level only the scalar field and not the metric 
fluctuation. It is also worthwhile to emphasize that the 
vacuum state defined by Sk^’ 7 * * /!)) = 0 is 100 % transla¬ 
tional and rotationally invariant (the formal proof was 
presented in Appendix A of Ref. [ 8 ]). 

Our next task is to connect the quantum theory of 
the inflaton perturbations with the primordial curvature 
perturbation. We proceed by choosing to work in the 
longitudinal gauge, and express Eq. 0, in terms of the 
expectation value of the conjugated momentum, this is, 

( 15 ) 


1 The Hankel functions of the first kind are defined as ll [) 1 (x) = 
J[/ (x) + iY v (x) with Ju and Y v the Bessel functions of the first 
and second kind respectively. 


It is clear that, in the vacuum state, {nu.{y)) = 0, which 
implies = 0 , i.e. there are no perturbations of the 
symmetric background space-time. It is only after the 
collapse has taken place (| 0 ) / | 0 )) that (dk (y)) e /0 
generically and $k / 0 ; thus, the primordial inhomo¬ 
geneities and anisotropies are born from the quantum 
collapse. 

It is also important to note that the quantum collapse 
affects all modes k of the inflaton, this is, the collapse 
takes the original vacuum state | 0 ) to a new quantum 
state: 

!©) = ••■ |B-k 2 )®|B-k 1 )®|0k o )<S>l©k 1 )®|0k 2 ) ■ ■ ■ (16) 

Given Eq. USD, which was provided by the semiclas- 
sical framework, and also that all modes of the inflaton 
field are now in a post-collapse state | 0 ), we can clearly 
see that the expectation value {^(y)) serves as a source 
for d/k for all k. Once the collapse has created all modes 
Tk, we can divide them in two types: 

1. Modes with an associated proper wavelength big¬ 
ger than the Hubble radius at the time of col¬ 
lapse, we will call these the super-horizon modes; 
i.e. their corresponding Fourier (conroving) modes 
satisfy k -C 

2. Modes with an associated proper wavelength 
smaller than the Hubble radius at the time of col¬ 
lapse, we will call these the sub-horizon modes; i.e. 
their corresponding Fourier (conroving) modes sat¬ 
isfy k^>'H{ql). 

In order to continue with the model, we will consider 
that the collapse is somehow analogous to an imprecise 
measurement of the operators y^' 1 (y) and tt r ’ 7 (y). This 
is, we need to specify the dynamics of the expectation 
values {y R,I {r))) and ( 7 T k ’ 7 (y )), evaluated in the post¬ 
collapse state. Furthermore, the analytical expression 
for {y R,I (v)) and {n RJ (y)) should depend on {y^ J (yl)) 
and ( 7 r k ’ J (^)}, i.e. the expectation values evaluated at 
the time of collapse. In the next subsection we will show 
the precise manner to specify (Ok’ 7 (Vk)) and R 

B. Collapse schemes 

At this point in the analysis we need to characterize in 
a more precise manner the collapse of the wave function. 


2 An imprecise measurement of an observable is one in which one 

does not end with an exact eigenstate of that observable but 

rather with a state that is only peaked around the eigenvalue. 

Thus, we could consider measuring a certain particle’s position 

and momentum so as to end up with a state that is a wave packet 

with both position and momentum defined to a limited extent 

and, which certainly, does not entail a conflict with Heisenberg’s 

uncertainty bound. 





Evidently, it would be desirable to provide a physical 
mechanism for the collapse. Nevertheless, there are some 
aspects that need to be addressed first. A full workable 
relativistic collapse mechanism is still unknown, however, 
some relativistic models have been recently proposed [26l . 
0 and are still under development. 

On the other hand, some non-relativistic objective col¬ 
lapse rnodels have been analyzed previously in the litera¬ 
ture Hunt- In particular, the Continuous Spontaneous 
Localization (CSL) model 0 is based on a non-linear 
stochastic modification of the standard Schrodinger equa¬ 
tion, in this way, spontaneous and random collapses of 
the wave function occur all the time, to all particles, re¬ 
gardless they are isolated or interacting. The idea behind 
the CSL model, sometimes referred to as the “amplifica¬ 
tion mechanism,” is that the collapses must be rare for 
microscopic systems, in order to not alter their quantum 
behavior as described by the Schrodinger equation, but 
at the same time, their effect must increase when several 
particles are hold together forming a macroscopic system. 

The CSL model has bee n ap plied to the inflationary 
universe in previous works flOl. I33l-I35l I45l| ; nevertheless, 
the results obtained in those works are different among 
each other (e.g. in Ref. 0 the amplitude of primordial 
gravitational waves is zero at first-order in the perturba¬ 
tions, meanwhile in Ref. (35l.l45j| this amplitude is similar 
to the predicted by the traditional inflationary model). 
The reason is the conceptual approach taken to address 
the subject, specifically, the treatment of the metric per¬ 
turbations (in 0 the metric is always classical, while 
in [35], HH metric perturbations are quantized). Further¬ 
more, there are still a few limitations on the CSL infla¬ 
tionary model that need to be investigated in detail, for 
instance: (i) CSL model is actually non-relativistic, but 
llOl I33H351 l45j assume a field theoretic CSL-like version, 
mode by mode, in momentum space, without a physi¬ 
cal justification; (ii) the CSL amplification mechanism 
is absent, or introduced as an ad hoc assumption. It is 
clear that further research is needed in order to consider 
a complete collapse mechanism and its successful imple¬ 
mentation to the inflationary universe. 

In the present manuscript, we are interested in ana¬ 
lyzing the characteristics of the observational predictions 
when considering a generic self-induced collapse. Thus, 
we will not consider a particular collapse mechanism and 
instead proceed in a pragmatical way. We will assume 
that whatever the collapse mechanism is behind, at the 
end of the collapse process, which can be associated to the 
time of collapse, we can characterize the expectation val¬ 
ues of the field and the momentum, evaluated at the post¬ 
collapse state. More precisely, we assume that the effect 
of the collapse on a state is analogous to some sort of 
approximate measurement; in other words, after the col¬ 
lapse, the expectation values of the field and momentum 
operators in each mode will be related to the uncertain¬ 
ties of the initial state. In the vacuum state, y k and ffk 
individually are distributed according to Gaussian wave 
functions centered at 0 with spread (Aj/k)o and (A#k)o> 


respectively. We could consider various possibilities for 
such relations; we will refer as “collapse schemes” to the 
different ways of characterizing the expectation values. 
In past works Ill three different schemes were consid¬ 
ered. These schemes were called independent, Newtonian 
and Wigner collapse schemes. In the following, we will 
describe them briefly. 


1. Independent scheme 

In this scheme one assumes that the expectation values 
''Ft I 

of the field’s mode y k ’ , and their conjugate momentum 
n^' 1 , acquire independent values randomly. The expec¬ 
tation at the time of collapse is assumed to be of the 
form 


(V k J ( ? 7k)) = \J ( A i'k’ / ( ? ?fe)) 0 > ( 17a ) 

<*k ’ J K)) = \j K)) 0 - ( 17b ) 

In this scheme the expectation value jumps to a ran- 
(Ft F) 

dom value ' multiplied by the uncertainty of the 
vacuum state of the field. The random variables j , 

are selected from a Gaussian distribution centered 
at zero, with unity spread, and are statistically uncorre¬ 
lated, that is the rationale of the name. This means that 
we are ignoring the natural correlation that exists in the 
conjugate fields in the pre-collapse state. InlBTl is shown 

the explicit form of (^Ay^’ 1 (r}^ and (^Att£ j (rify'j 
within this collapse scheme. 


2. Newtonian collapse scheme 

This scheme is motivated by the fact that in the equa¬ 
tion for the Newtonian potential, Eq. only the ex¬ 

pectation value of 7Tk appears. Thus, one is led to con¬ 
sider a scheme where the collapse affects only the conju¬ 
gated momentum variable, this is 


<2/k J (^)> = 0, 



(18) 

As in the previous case, x'^, 2 r> represents a random Gaus¬ 
sian variable normalized and centered at zero. The quan¬ 
tity (^Air^’ 1 within this collapse scheme is the same 

as in the independent scheme (see IBID . 
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3. Wigner collapse scheme 

The last collapse scheme, analyzed in detail in Refs. 
I i, is motivated by considering the correlation be¬ 
tween y R ' 1 and n R ’ 1 existing in the pre-collapse state 
and characterize it in terms of the Wigner function (one 
knows from Heisenberg’s uncertainty principle that the 
field and momentum variables should be correlated). 

The Wigner function of the vacuum state is a bi- 
dimensional Gaussian function. The assumption is that, 
at a certain (conformal) time 77 ^, the part of the state 
characterizing the mode k will collapse, leading to a new 
state in which the expectation value of the fields will be 
characterized by 

(Vk 1 ^)) = x k J AkiVk) cos&k{Vk), (19a) 

<^k J (Vk)) = x k 1 kh-k(rik ) sin ®fc( ? ?fc)> (19b) 

where x^’ 1 is a random variable, characterized by a 
Gaussian probability distribution function centered at 
zero with spread one. The parameter A k(Vk) represents 
the major semi-axis of the ellipse characterizing the bi- 
dimensional Wigner function that can be considered a 
Gaussian in two dimensions; this is, the ellipse corre¬ 
sponds to the boundary of the region in “phase space” 
where the Wigner function has a magnitude larger than 
1/2 its maximum value. The other parameter, repre¬ 
sented by 0 fe(? 7 fe), is the angle between that axis and the 
y k ’ axis. The details involving the Wigner function and 
the collapse scheme can be consulted in Ref. Q- The 
parameters A k and 0 fc depend on the time of collapse, 
one can follow the analysis presented in Ref H] in order 
to find an expression for A k and 0 & in terms of the time 
of collapse, but bearing in mind that, for the present 
manuscript, we are considering a quasi-de Sitter universe 
[see EH for an explicit expression of A k(Vk) an d 0 fc(? 7 fc)] 


C. The curvature perturbation for the three 
collapse schemes 

So far, we have established the relation between the 
curvature perturbation and the quantum matter fields 
[see Eq. (ITol) ] and characterized the collapse by means 
of the expectation values of the field and its momentum, 
i.e. by introducing the collapse schemes. The next aim 
is to present an explicit expression for the curvature per¬ 
turbation in terms of the parameters characterizing each 
collapse scheme. In order to attain that goal, we must 
first find an expression for the evolution of the expecta¬ 
tion values of the fields. In fact, as can be seen from 
Eq. (HU), we will only be concerned with the expectation 
value of the conjugated momentum (7Tk(?7)). In IB 21 we 
show that 


= F{kr],z k ){y* J (ril))e 
+ G(A;j7,^)(w k ’ / (^)) e , (20) 

with the definitions of F(ky,Zk ) and G(krj,z k ) also in 
IB 21 The parameter z k is defined as Zk = ki 7 ^; thus, Zk is 
directly associated to the time of collapse 77 ^. 

Finally, substituting Eq. m in Eq. ms , we can 
find an expression for the curvature perturbation (in the 
longitudinal gauge). 


^k(r?) 


f^v H 

V 2 Mpk 2 


X 


F(h 7, z k ) ({y£(vk))e + i(yL(Vk))e) 


+ G{kr 7 , Zk) ((fr k (?7fc))e + *(dk(f?fc)}e) 


■ ( 21 ) 


We can see from Eq. m how the curvature pertur¬ 
bation depends on the three collapse schemes through 
the quantities (2/k( r 7fc))e and (^k^/-))© • Henceforth, we 
have three different expressions for *Fk(? 7 ) corresponding 
to the three previously introduced collapse schemes char¬ 
acterized in Eqs. (fT7l) . (fl8ll . (fl9l) (explicit expressions of 
d^k in the three schemes are given in IB 21) . 

One useful gauge-invariant quantity often encountered 
in the literature is the variable TZ(x). The field TZ(x) 
is a field representing the curvature perturbation in the 
comoving gauge. Its Fourier’s transform, represented by 
IZk, is constant for modes “outside the horizon” (irre¬ 
spectively of the cosmological epoch), i.e. for modes with 
k <C H = aH (and assuming adiabatic perturbations). 
This is, the value of TZk during inflation (in the limit 
k <C ki) would remain unchanged at all times, until the 
mode “re-enters the horizon,” namely when k ~ 'H. 

On the other hand, the curvature perturbation *F in 
the longitudinal gauge, is also constant for modes outside 
the horizon during any given cosmological epoch but not 
during the transition between epochs. In fact, during 
the transition from the inflationary stage to the radia¬ 
tion dominated stage, 'F is amplified by a factor of l/ey 
[iiL [46j. This behavior differs from the one of TZ, which 
remains constant in spite of the epoch transition. 

The curvature perturbation in the comoving gauge TZ 
and the curvature perturbation in the longitudinal gauge 
^F are related as 7 Z = ^F + (2/3)('H^ lv F' + \F)/(1 + w), 
with uj = p/p. Therefore, for modes such that k <C kL, 
during the inflationary epoch w + 1 ~ 2ey/3, one has 


lim i?k — lim 


ev 


( 22 ) 


with ^( 77 ), calculated during inflation, in the limit such 
that the modes are well outside the “horizon” (i.e. in the 
regime where |fc? 7 | <C 1). 
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Therefore by expanding the expressions 'I'k, within the 
three collapse schemes, to the lowest order in \krj\, and 
by making use of Eq. (l22l) . we can find 77k- Thus, af¬ 
ter performing such expansion, the comoving curvature 
perturbation is 


77^ ~ R k 


M(\z k \)X K i+N(\z k \)X K2 


\kr/\ 


3/2-v 


(23a) 


77Jf wt ~ R k N(\z k \)X K2 \kr 1 \ 3 ' 2 - 1 ', (23b) 


P(fc) = A s fc ns_1 , where A s is the amplitude and n s is 
known as the spectral index of the scalar perturbations. 
Thus, in this section we will find a similar expression in 
which we will identify the amplitude and the scalar index 
within the collapse model. 

We begin by showing how the observational quantities 
can be related with the parameters characterizing the col¬ 
lapse. The temperature anisotropies ST/Tq of the CMB 
are clearly the most direct observational quantity avail¬ 
able, with To the mean temperature. Expanding ST/Tq 
using spherical harmonics, the coefficients ai m are 


77™ g ~ R k W{\z k \)Xk\kr}\ 3/2 - v . (23c) 

The functions M(| Zk |), 7V( \Zk\) and W(\Zk \) are defined 
in IB 31 and the amplitude 
Rk = sjL^lt v H2 v ~ 11 l' 1 r(v - 1 )/M P fc 3 * * * / 2 ; 
also we have introduced the definitions Zk = kr ]£ and 
A k = a-'k + ix k- 

Equations (l23l) are the main result of this section. 
They relate the initial curvature perturbation, which 
is associated with the temperature anisotropies in the 
CMB, with the parameters characterizing each collapse 
scheme, i.e. the time of collapse and the random vari¬ 
ables. There is no analogous expression in the tradi¬ 
tional inflationary paradigm, in which by relying on some 
“quantum-to-classical” arguments (see Refs, [lj, [2j] for a 
detailed discussion on the conceptual problems regard¬ 
ing such arguments), one is able to go from 77k to 77k 
but without a clear identification of the physical (and 
probably random) process that originated the classical 
curvature perturbation. 

We strongly remark that the random variables corre¬ 
sponding to each collapse scheme are fixed after the col¬ 
lapse of the wave function has occurred. In other words, 
if we somehow knew their exact value, we would be able 
to predict the exact value for 77k; notice that we have not 
even mentioned notions such as average over an ensemble 
of universes or some related concepts. Nevertheless, we 
will do make use of the statistical properties of the ran¬ 
dom variables to be able to make theoretical predictions 
for the observational quantities, e.g. the power spectrum 
and the spectral index; this will be the focus of the next 
section. 


IV. AN EQUIVALENT POWER-SPECTRUM 
FOR THE CURVATURE PERTURBATION 

The focus of this section is to find an equivalent expres¬ 
sion for what it is commonly referred to as the primordial 
power spectrum for the scalar perturbations. In the stan¬ 
dard inflationary paradigm, such expression is given bj0 


3 Actually, in the literature, one finds two kinds of power spec¬ 

trum: the dimensional power spectrum ’P(fc) and the dimension¬ 

less power spectrum P(fe); the latter is defined in terms of the 

former by P{k) = (k 3 /2ir 2 )V(k) 


aim = J €>(h)Y l * m (0,ip)dQ,, (24) 

with n = (sin 6 sin ip, sin 9 cos tp, cos 9) and 9 , tp the coor¬ 
dinates on the celestial two-sphere; we have also defined 
0(n) = 5T(fi)/T 0 . Assuming instantaneous recombina¬ 
tion, the relation between the primordial perturbations 
and the observed CMB temperature anisotropies is 

1 f Va 

0(n) = [T + -5 7 ](??D)+n-fT 7 (? ?D ) + 2 / ^'(r?)dry, (25) 

4 Jri D 

where tjd is the time of decoupling; 5 7 and Vj are the 
density perturbations and velocity of the radiation fluid. 

It is common practice to decompose the temperature 
anisotropies in Fourier modes 

0 (^) = E (26) 

k 

with Rp the radius of the last scattering surface. Af¬ 
terwards, one solves the fluid motion equations with the 
initial condition provided by the curvature perturbation 
during inflation. Furthermore, using that e lk ' Rnn = 
47r Xwm i l ji{kRD)Yim{9,p)Y ] k m (k) 1 expression (O can be 
rewritten as 

aim = -jyY.j^ kR D)YL{km k), (27) 

k 

with ji(kRp) the spherical Bessel function of order l. 

The linear evolution that relates the initial curvature 
perturbation IZk and the temperature anisotropies 0(k) 
is summarized in the transfer function T(k), in other 
words, T(k) is the result of solving the fluid motion equa¬ 
tions (for each mode) with the initial condition provided 
by the curvature perturbation 77-k and then make use of 
Eq. (l25l) to relate it with the temperature anisotropies. 
Thus 0(k) = T{k)TZ k . 

Consequently, the coefficients a; m , in terms of the 
modes 77k, are given by 

aim = ^-J2^ kR D) Y im(k)T(k)n^ (28) 

k 
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with 7?-k during inflation, and in the limit k <C R or 
equivalently \kr ]| -C 1. 

By substituting Eqs. E51) . corresponding to the ex¬ 
plicit form of 7?-k for each collapse scheme, in Eq. (j28| , 
one can see how the coefficients ai m are directly related to 
the random variables X k. Consequently, the coefficients 
aim are in effect a sum of random complex numbers (i.e. 
a sum over k where each term is characterized by the 
random variables X k, which is itself a complex random 
number), leading to what can be considered effectively as 
a two-dimensional (i.e. a complex plane) random walk. 
As is well known, one cannot give a perfect estimate for 
the direction of the final displacement resulting from the 
random walk. Nevertheless, one might give an estimate 
for the length of the displacement for the random walk. 
Such length is naturally associated with the magnitude 
\aim\ 2 ', hence, we can make an estimate for the most likely 
value of \ai m \ 2 and interpret it as the theoretical predic¬ 
tion for the observed \ai m \ 2 . Moreover, since the collapse 
is being modeled by a random process, we can consider 
a set of possible realizations of such process characteriz¬ 
ing the universe in an unique manner, i.e., characterized 
by the random variables X^. If the probability distribu¬ 
tion function of X^ is Gaussian, then we can identify the 
most likely value |azm|ML w ith the mean value \ai m \ 2 of 
all possible realizations; this is, = |a; m | 2 . The 

most likely value |az m ||j L in each collapse scheme is ex¬ 
plicitly given in IHl 

Moreover, we need to make some further assumptions 
regarding the ensemble average of the product of the ran¬ 
dom variables for each collapse scheme. We will assume 
that x k ’ variables are uncorrelated, this is, they satisfy 


x k s X k' s = ^k.k' + ^k,-k'; ^k s X k' s = ^k.k' — ^k,-k'; 

.. .. ’ (29) 

with the label s indicating the particular scheme associ¬ 
ated to the random variables. 

Note that we have taken into account that the real and 
imaginary parts of the random variables are independent 
for every scheme. Additionally, we have considered the 
correlation between the modes k and k in accordance 
with the commutation relation given by [a k , a k 7] and 

[«k»«k'] I see Ec h ©]• 

The observational data is presented in terms of a quan¬ 
tity Ci called the angular power spectrum. The def¬ 
inition of Ci is given in terms of the coefficients ai m 
as Ci = (2/ + l) -1 |«Zm| 2 - Therefore, we can use 

the prediction for |azm|ML fo r each collapse scheme and 
give a theoretical prediction for Ci for the three collapse 
schemes considered. Thus, Eq. (l29l) and using |az m |ML 
for each collapse scheme, we obtain 

r°° rik r 

Ci=4n —jf(kR D )T(k) 2 —Q(\zk\)k 3 ~ 2u , (30) 

Jo & n 

where 


C = (2 v ~ n/2 T(v - l)H M 3 / 2 -") 2 , (31) 

and we have taken the limit L —» oo and k — > continuum 
in order to go from sums over discrete k to integrals over 
k. The function Q{\z k \) varies for each collapse scheme. 
For the independent scheme 


Q(|* fe |) ind = M 2 (|^|)+iV 2 (|^|); 

(32) 

for the Newtonian scheme 


Q(\z k \) newt = N 2 {\z k \), 

(33) 

and for the Wigner scheme 


Q(N) wig = W 2 (|* fc |); 

(34) 


the definitions of M, N, W are in IB 31 

In the standard inflationary paradigm, a well-known 
result is that the power spectrum P(k) for the perturba¬ 
tion 7?-k and the Ci are related by 

r °° ji. 

Ch = to] —jf(kR D )T(k) 2 P(k). (35) 

Thus, by comparing Eq. (l30l) with (f35l) we can ex¬ 
tract an “equivalent power spectrum” for each collapse 
scheme (more details can be found in iDl) . The form of 
the power spectrum, within the three collapse schemes, 
has a generic form 

P{k) = 4 Q(\zk\)k 3 ~ 2v , (36) 

7T Z 

Equation (l36ll is the main result of this section. In 
the next section, we will compare our prediction with 
the one given by the traditional approach and make a 
couple of remarks regarding the physical implications of 
our prediction. 

We would like to end this section by noting that our 
prediction for the power spectrum was extracted from 
what in principle are observable quantities, i.e the Ci s. 
In fact, our model gives a direct theoretical prediction 
for the observed Ci, Eq. (E}, and then from such ex¬ 
pression we “read” what can be identified as the “power 
spectrum” in the traditional approach of inflation. This 
is conceptually different from the orthodox approach [2lJ 
in which the power spectrum is obtained directly from 
the two-point correlation function (0|'k(k)'&(k , )|0). In 
contrast, our power spectrum is obtained from (irkKifk'), 
where the expectation values are evaluated at the post¬ 
collapse state, this is, we have never relied on the calcu¬ 
lation of the two-point quantum correlation function. In 
m we show in detail the calculation of the CMB temper¬ 
ature angular power spectrum and its relation with the 
scalar power spectrum, this serves to present the reader 
why our proposal does not rely on the quantum two-point 
correlation function. 
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V. COMPARISONS BETWEEN THE 
TRADITIONAL AND THE COLLAPSE POWER 
SPECTRUM 


Let us recall that the standard prediction for the scalar 
power spectrum for the curvature perturbations, within 
single field slow-roll inflation, is P(k) = Ag/c" 8-1 , where 


2 2p - 4 \T{u)\ 2 H 2 \r]\ 3 - 2i> 

ir 3 MpCH 


n s — 1 = —6 en + 2 Sy, 


. (37) 

with v = 3/2 + Sen — Sy. On the other hand, within our 
model, v = 3/2 + en — Sy [see Eq. (fl2l) ]: therefore, the 
equivalent scalar spectral index is — 1 = —2e# + 2 Sy. 

As has been analyzed in Ref. |47|, the quantity A s , 
namely the amplitude, is approximately a time indepen¬ 
dent function [i.e., d/di^{H 2 \ri\ 3 ~ 2 ' y /en} = OiejjAy)] for 
all 77, this is, even if H , eu and \r]\ 3 ~ 21 ' are time depen¬ 
dent, their combination, as it appears in A s , is essen¬ 
tially constant. Since A s is nearly time independent, it 
is customary to express P(k) in terms of the value of the 
conserved quantities when the mode crossed the horizon, 
k = aH. In other words, one chooses to express the value 
of the power spectrum, which is written in terms of a con¬ 
served quantity at first order in the slow-roll parameters, 
as a time independent function of wavenumber k , 


where 


P(k) 


2 2 ^~ 4 |r(i>)i 2 Hl(k) 

7 r 3 M 2 e*.(Jfe) 


(38) 


wavelength A p = a/k (with k the comoving wavenum¬ 
ber), at the time of collapse, is smaller or bigger than 
the Hubble radius H~ 3 , which we know is approximately 
constant during inflation. The answer is simply to com¬ 
pare which of the inequalities gets satisfied k 
or k -C a(r]/.)H, if it is the former then the mode is still 
“inside the horizon,” and we know that it will eventually 
cross the horizon, during inflation, and then “freezes;” 
however, if it is the latter then the mode is “bigger than 
the horizon” and is already “frozen,” consequently for 
this last type of modes it would make no sense to evalu¬ 
ate the power spectrum at the “horizon crossing.” 

The second important aspect of the collapse power 
spectrum is the function Q(\zk\), which is in principle a 
function of k. If one assumes that the time of collapse is 
of the form 77 J] oc fc , then Zk = kr = z is independent 
of fc; consequently, < 3 (|z|) would be also k independent; 
thus, the collapse power spectrum would depend on k as 
P{k) oc k 3 - 2 ", which for all phenomenological purposes 
would be indistinguishable from the prediction given by 
the traditional inflationary approach. Consequently, in 
our model, if the dependence on k of the time of collapse 
77 )) is slightly different from 77 £ oc 1/fc, then our proposal 
will yield different predictions from the standard infla¬ 
tionary paradigm. 

Summarizing the above discussion, if 77 ^ oc A; , the 
collapse power spectrum can be expressed as P(k) = 
Ak Us ~ 1 : with the amplitude given by 








H 2 (k) ( H 2 \ V \ 3 ~ 2 > \ 

eM fc ) V £h ) 


k—aH 


(39) 


then, one computes the spectral index as n s — 1 = 
d\nP(k)/d\nk using Eqs. (l38l) and (l39l) . 

On the other hand, in the collapse model, the equiva¬ 
lent power spectrum is given by P(k) = C /rr 2 Q{\zk\) 
k 3 ~ 2u , Eq (l36ll . A few remarks are in order: 

First, the quantity C is of the same structure as A s 
111 the traditional approach [see Eq. ( 1 ■ 1 1 1) and Eq. ©]■ 
In other words, is nearly a time independent function, 
that is to say its derivative with respect to the confor¬ 
mal time is of second order in the slow-roll parameters. 
One could follow the standard procedure and re-express 
C when the mode crossed the horizon; however, within 
our approach, the curvature perturbation and, thus, the 
power spectrum is non-vanishing only after the time of 
collapse. 

Specifically, once the collapse has occurred (or more 
precisely, the collapse mechanism has ended and the 
semiclassical approximation is valid) at some time 77 J] 
during inflation, and as a consequence, Tk 7 ^ 0 for 
all k, i.e. the primordial curvature perturbation has 
been generated. Then, one can focus on some partic¬ 
ular Fourier’s mode Tk and ask if its associated proper 


we emphasize, that A is constant in time, up to second 
order in the slow-roll parameters, and independent of fc; 
additionally, the scalar spectral index obtained from the 
collapse model is 


n.- 1 = -2e ff + 2S V , (41) 

which is a little different from the standard prediction 
[see Eq. (1371) ]. 

One can relax the condition on the form of the time 
of collapse by allowing a small dependence on k, namely, 
assuming that the time of collapse is 77 ^ = A/k + B and 
use the observational data to constraint the parameters 
A and B (a similar analysis has been done for the perfect 
de Sitter universe, see Ref. 0 )- In such case, the collapse 
power spectrum is of the form P(k) = C/'rr 2 Q(\zk\)k 3 ~ 2 '', 
with Zk = A + Bk [Eq. (1M1) ]. This is, the collapse power 
spectrum depends on k in two ways, one as k 3 ~ 2l/ and 
another through the time of collapse Zk = kr = A + Bk; 
however, the former dependence is due to the inflation¬ 
ary dynamics of the curvature perturbation, and, thus, 
interpreted as the scalar spectral index, while the lat¬ 
ter reflects the consideration of the collapse hypothesis. 
Additionally, one could also reinterpret the dependence 
on k, introduced by the collapse proposal through the B 
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parameter, as something resembling to a running of the 
scalar spectral index; this would be an interesting idea 
to pursuit since it has been pointed out before that the 
lowest multipoles of the temperature anisotropies prefer 
models with n lun 7 ^ 0 0|. However, the physical 

interpretation of the departure of a power law in k corre¬ 
sponding to the primordial power spectrum is different. 
In the collapse proposal, it arises from the inclusion of 
the self-induced collapse of the wave function while in 
the standard approach it arises from considering the sec¬ 
ond order approximation in the slow roll parameters (see 

Ref.IH) 

As mentioned before, the time of collapse can occur at 
any time during the inflationary regime; in particular, it 
can happen when the proper wavelength of the mode is 
bigger or smaller than the Hubble radius, which is ap¬ 
proximately constant. 

If the proper wavelength of the mode, at the time of 
collapse, is bigger than the Hubble radius, then k <C 
a{i^)H = H^) cs —1 /r 7 fc which is equivalent to — kr] £ 

•C 1. Then, for modes with a proper wavelength big¬ 
ger than the Hubble radius at the time of collapse, the 
collapse power spectrum can be approximated by 

P(k) ~ ^E(\z k \)k 3 - 2 ", (42) 

TT 2 

where S(| 2 fc|) is obtained by expanding the function 
Q(\z k \), when \z k \ -A 0 [i.e., expanding M(\z k \),N(\z k \), 
and W(|zfc|)], to the lowest order in \z k \. Thus, for each 
scheme the function E(\z k \) is 


s(N) ind ^ 


4 





(43a) 


N) newt = 


Zk 


4tt 2 u 2 (p — l ) 2 


(43b) 


;(M) wig = 


16 


\z k 


1 

4|zfc| 

1/-1/2 


\ + ( x ~ Vi + 10N 2 + 9|**| 4 ) 

1,+1/2 


-1 2 


[r(v)2 


1/-1 


COS 01 


Zk 


T{v — 1)2' 


■ sin 0 i. 


(43c) 


where tan20fc ~ 4\z k \/(l — 3|zfc| 2 ). 

On the other hand, if the proper wavelength associated 
to the mode is smaller than the Hubble radius, at the 
time of collapse, then k a(r/^.)H, which is equivalent 
to — k'q’f. /§> 1. Then the approximated collapse power 
spectrum, when —kr)%. = \z k \ —> 00 , is 


P{k) =* ^T(\z k \)k 3 - 21 ', (44) 


where T(|^|) is obtained by considering the asymptotic 
behavior of the function < 5 (|zfe|) [i.e., the asymptotic 
behavior of M(\z k \),N(\z k \),W(\z k \)] when \z k \ -A 00 . 
Thus, for each scheme the function T(|^fc|) is 


T(|~fc|) ind = - 

n 2 


1 + 


1 


z k \ 

1 2 


sin fi(v, \z k \) + 

r(i/ + 5/2) \ 

2T(j/+1/2)7 J 


x( 2r , r(^ + 5/2) \ 
1 + 2T(i/+1/2)J 


r(i/ + 3/2) V 
4 |^ fe | 2 \T(v- 1 / 2)7 

cos/3(u, \z k \) 


1 + 


1 




cos/3(r+ | z k j) - 


sin/3(i/, |z fc |) T(v + 3/2) 
2\z k \ r(v- 1 / 2 ) _ 


,(45a) 


T(|+fc|) newt = 


1 + 


\z k 


7T 

r(i + V2)_V 

2 ' 2r(i/ + 1/2 )) 


cos /3(v, \z k \) - 


sin |z fc |) T(v + 3/2) ' 
2\z k \ T(v - 1/2) 


,(45b) 


T(|z fe |) wig = ^ 

7 T 2 


2v 


\z k 


13/2 


x ^cos f3(v, \z k \) - 
- ^sin/3(i/, \z k \) + 
+ cosf3(v, \z k \) — 


2 \z k 

!>- 1/2) 

cos /3(zx, 

z k \)T(v + 5/2) 

2 | z k 

r^ + 1/2) 

smP(u, | 

z k \)T(u + 3/2)1 

2 | z k \ 

r(zx-i/ 2 )J 


COS 0fc 
2 


sin 0 ^ 


(45c) 


where f3(v, \z k \) = \z k \ — (-k/2){v + 1/2) and tan20£ ~ 
-4/3|**|. 

Expressions m and (1451) are useful for performing 
the comparisons between the theoretical prediction of our 
model and the observational data. 


VI. PRIMORDIAL POWER SPECTRUM IN 
QUASI-DE SITTER 

In this section, we will show the primordial power spec¬ 
trum in the quasi-de Sitter case for the different collapse 
schemes analyzed in this paper. 

In particular, we will focus on the cases in which the 
proper wavelength associated to the mode is bigger and 
smaller than the Hubble radius at the time of collapse, i.e 
we analyze the cases such that the comoving k , associated 
to the curvature perturbation, is bigger or smaller than 
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a(r]£)H = 7i(j]1)\ these cases correspond to Eqs. (l42l) 
and (TH1) respectively. 

This preliminary qualitative analysis indicates that the 
aforementioned collapse schemes are good candidates to 
account for the observational data of the CMB collected 
by the Planck [28| and WMAP [29| collaborations. How¬ 
ever, we will not perform here the statistical analyses 
to compare the theoretical predictions with the observa¬ 
tional data in order to constrain the value of the free 
parameters of the collapse model ( A and B). We will 
leave this task for a forthcoming paper jjl|. 

First, let us define a fiducial model with a primordial 
power spectrum P(k ) = A s k ns ~ 1 with n s = 0.96 that 
will be taken just as a reference to discuss the plots we 
obtain for the collapse models. The value of n s , for our 
fiducial model, is the mean value obtained by the Planck 
and WMAP collaborations. Let us remind the reader 
that the free parameters of all collapse schemes (A and 
B) are related to the time of collapse of each mode of the 
inflaton field = A/k + B. 

Figures in m and [3] show the primordial spectra for 
the different collapse schemes ( independent, Newtonian, 
Wigner ), in the case where the proper wavelength asso¬ 
ciated to the mode is bigger than the Hubble radius at 
the time of collapse, i.e. k <C B(r^) [Eqs. ([43]) ]. for dif¬ 
ferent values of the collapse parameters A and B. The 
primordial power spectrum of the fiducial model is also 
shown in each Figure; the value of A s (detailed in the 
caption of each Fig.) is settled in each case in order 
to provide the reader a clear idea regarding the differ¬ 
ences in the form of the different spectra shown in the 
plot. For example, in Fig. [1] the fiducial power spectrum 
could be normalized in such a way that it almost over¬ 
laps the primordial spectrum for B = —0.1 Mpc. How¬ 
ever, for B = —0.5 Mpc and B = — 1 Mpc there is no 
value of A s that makes both spectra (collapse and fidu¬ 
cial model) overlap. Hereinafter, we discuss the effects 
of introducing the collapse hypothesis for the different 
collapse schemes in the primordial power spectrum. The 
relevant values for k that will affect the prediction of 
the CMB temperature and polarization anisotropy are 
10 -6 Mpc -1 < k < 10 -1 Mpc -1 . Therefore, we have 
drawn a vertical line in each Figure, at k = 0.1 Mpc -1 . 

Figure [T] shows no change in the slope of the primordial 
power spectrum (with respect to the fiducial model one) 
for the independent scheme in the range k < 0.1 Mpc -1 
and the considered values of B and the two values of A; a 
notorious deviation in shape from the fiducial one occurs 
for k > 0.1 Mpc -1 . However, we have already men¬ 
tioned that the relevant values of l that affect the CMB 
spectrum correspond to 10 -6 Mpc -1 < k < 10 -1 Mpc -1 , 
thus we expect no deviation of the CMB temperature 
spectrum respect the fiducial one for this scheme. 

Again, Figure [2] shows no change in the slope of the 
power spectrum, but unlike the previous case, there is a 
deviation from the fiducial power spectrum for increas¬ 
ing values of B , and this deviation happens within the 
relevant range of k. Furthermore, it should be noted the 




Figure 1. Primordial power spectra for the independent 
scheme in the case where k <C T-L{ 77 ^). The power spectrum of 
the fiducial model is also shown in red. Different values of the 
collapse time = A/k + B are considered, the scalar spectral 
index n s = 0.96; Top: A = —10 1 , Bottom: A = —10 -2 ; for 
both figures A s = 0.39. The difference in the amplitude be¬ 
tween the collapse models and the fiducial model is artificially 
set to show the functional form of both models. 


huge value of A s considered in order to make compar¬ 
isons between the collapse and fiducial power spectrum. 
A different value of this parameter with respect to the 
the value of the standard model could result in a differ¬ 
ent value for the scale energy of the inflationary period. 
However, the complete determination of this scale energy 
depends not only on the temperature fluctuations of the 
CMB but also on the B-mode polarization. As men¬ 
tioned before, a recent joint paper by the BICEP2/Keck 
and Planck collaborations concludes that there is no evi¬ 
dence for primordial B-modes at low angular multipoles. 

Figure [3] shows the differences in the primordial spec¬ 
tra among different values of B for the two values of A 
considered and also with respect to the fiducial model 
for the Wigner’s scheme. It is interesting to note that 
the deviation with respect to the fiducial model goes to 
the opposite side that in the Newtonian scheme discussed 
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Figure 2. Primordial power spectra for the Newtonian col¬ 
lapse scheme in the case where k -C %(??(;). The power spec¬ 
trum of the fiducial model is also shown in red. Different 
values of the collapse time rfe = A/k + B are considered, 
the scalar spectral index n s = 0.96; Top: A = — IGF 1 and 
A a = 2.8 x 10 s , Bottom: A = —1CF 3 and A g = 3 X 10 13 . The 
difference in the amplitude between the collapse models and 
the fiducial model is artificially set to show the shape of both 
models. 


previously. Furthermore, as in the Newtonian scheme, it 
was necessary to consider large values of A s in order to 
make the collapse and fiducial model spectra comparable. 
Thus, the discussion made for the Newtonian scheme is 
also valid for this case. 

In summary, in this preliminary analysis for the case 
where k <C we could find values of the collapse 

parameters ( A and B) that make the primordial power 
spectrum of the collapse models to be almost equal to 
the fiducial primordial power spectrum for the all the 
schemes considered in this paper. In consequence, we 
expect that the corresponding CMB spectrum will not 
deviate much from the fiducial one. 

Figures |5J [5] and [G] show the primordial spectra for 
the different collapse schemes ( independent, Newtonian, 
Wigner ) in the case where the proper wavelength asso- 


Figure 3. Primordial power spectra for the Wigner collapse 
scheme in the case where k -C The power spectrum 

of the fiducial model is also shown in red. Different values 
of the collapse time rg. = A/k 4- B are considered, the scalar 
spectral index n s = 0.96; Top: A = —10 -1 with A a = 80, 
Bottom: A = —1CF 2 with A a = 7000. The difference in the 
amplitude between the collapse models and the fiducial model 
is artificially set to show the functional form of both models. 


dated to the mode is smaller than the Hubble radius at 
the time of collapse, i.e. k H(rff. ) [Eqs. (l45l) ] for dif¬ 
ferent values of the collapse parameters A and B. The 
primordial power spectrum of the fiducial model is also 
shown in each Figure and the value of A s (detailed in 
the caption of each Figure) is settled in each case in 
order to provide the reader a clear idea regarding the 
differences in the form of the different spectra shown 
in the plot. The relevant values for k that will affect 
the prediction of the CMB temperature and polarization 
anisotropy are 10 _6 Mpc _1 < k < 10 -1 Mpc _1 . There¬ 
fore, we have drawn a vertical in each Figure, at k = 0.1 
Mpc -1 . For the independent scheme (Fig. 0}, it should 
be noted that for A = —10 2 the power spectrum deviates 
from the fiducial model for increasing values of B , how¬ 
ever, it is a small change compared to the deviations of 
other schemes (see Figs. [5]and[n|), while for A = — 10 6 
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Figure 4. Primordial power spectra for the independent 
scheme in the case where k The power spectrum of 

the fiducial model is also shown in red. Different values of the 
collapse time = A/k + B are considered, the scalar spec¬ 
tral index n a = 0.96; Top: A = —10 2 , Bottom: A = —10 6 , for 
both figures A a = 0.43. The difference in the amplitude be¬ 
tween the collapse models and the fiducial model is artificially 
set to show the shape of both models. 


there is no difference in the spectrum among different val¬ 
ues of B. This is due to the fact that for large values of 
A 1 the value of zk becomes also very large, and therefore 
the value of B does not affect the final form of the power 
spectrum. This behavior is similar to the case where the 
scale factor is exactly de Sitter (see Ref. 0 )- 

For the Newtonian scheme (Figure[5]) it should be men¬ 
tioned that in both cases ( A = —10 2 and A = —10 5 ) 
considered for increasing values of B the spectrum devi¬ 
ates from the fiducial one for k > 0.01 Mpc -1 . For the 
Wigner scheme (Figure O, it should be mentioned that 
for increasing values of B the power spectrum deviates 
more from the fiducial model in both cases (A = —10 2 
and A = —10 6 ) for k > 0.005 Mpc -1 . In summary, in¬ 
creasing the value of B results in a deviation from the 
fiducial model spectrum; the deviation is more drasti¬ 
cally for the Newtonian and Wigner scheme. 




Figure 5. Primordial power spectra for the Newtonian col¬ 
lapse scheme in the case where k T-L{i 7 ^). The power spec¬ 
trum of the fiducial model is also shown in red. Different 
values of the collapse time rf/. = A/k + B are considered, 
the scalar spectral index n 3 = 0.96; Top: A = — 10 2 with 
A s = 0.3, Bottom: A = —10 s with A s = 0.42. The difference 
in the amplitude between the collapse models and the fiducial 
model is artificially set to show the shape of both models. 


VII. PREDICTIONS OF THE COLLAPSE 
SCHEMES ON THE CMB TEMPERATURE 
SPECTRUM 

The aim of this section is to show, that introducing the 
collapse of the inflaton wave function during inflation has 
observable consequences on the CMB fluctuation spec¬ 
trum. In this paper, we will limit ourselves to the analysis 
of the temperature auto-correlation spectrum; however, 
from a previous analysis of similar models 0 we might 
expect that the .E-mode polarization and temperature- 
E-mode cross correlation will be also modified as a con¬ 
sequence of the collapse hypothesis. As we will see, the 
effect is different for the three collapse schemes proposed 
in this paper and it also depends on the value of the 
time of collapse. This is not a surprise since we have 
shown in the previous section that there are differences 
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Figure 6. Primordial power spectra for the Wigner collapse 
scheme in the case where k W(r/£). The power spectrum 
of the fiducial model is also shown in red. Different values 
of the collapse time rff. = A/k + B are considered, the scalar 
spectral index n a = 0.96; Top: A = — 10 2 with A a = 0.39, 
Bottom: A = —10 6 with A a = 0.21. The difference in the 
amplitude between the collapse models and the fiducial model 
is artificially set to show the shape of both models. 


between the primordial power spectrum in each collapse 
scheme and the standard inflationary model one. We 
want to stress that, in this paper, we will only perform a 
preliminary analysis of the CMB spectrum predicted by 
the collapse models for some particular cases; a complete 
analysis, including statistical analysis, in which recent 
CMB data are confronted with the predictions from all 
collapse schemes, is in progress a In order to perform 
our analysis, let us define the cosmological parameters of 
our fiducial model: baryon density in units of the critical 
density fish 2 = 0.02212, dark matter density in units 
of the critical density ftcDMh 2 = 0.1187, Hubble con¬ 
stant in units of Mpc _1 kms _1 Hq = 67.75, reionization 
optical depth, r = 0.092, and the scalar spectral index, 
n s = 0.96. These are the best-fit values presented by 
the Planck collaboration [52 j. The value of A s is settled 
in each case in order to match the maximum of the first 


Doppler peak with the fiducial model one. 

Figs. [3 [8] and [9] show the prediction for the CMB 
temperature fluctuation spectrum for the independent, 
Newtonian and Wigner schemes in the case where k <C 
TL(r]^). Let us first analyze the case of the indepen¬ 
dent scheme; for the two cases shown in this paper 
(. A = —10 _1 and A = — 10 -2 ) there is only a very tiny dif¬ 
ference between the fiducial spectrum and the one includ¬ 
ing the self-induce collapse of the inflaton’s wave func¬ 
tion. We have calculated the x 2 using the latest release 
of Planck data, and the difference between the value of 
the fiducial model and the collapse model is within the 
expected values for this quantity. In consequence, we ex¬ 
pect that any value of B can explain recent observational 
data. 

Conversely, Figs. [5] and [9] show large deviations for 
increasing values of B in the temperature power spec¬ 
trum respect to the fiducial one. On the other hand, 
we would like to stress that the values of A s required to 
match the collapse and the fiducial spectrum are large 
and, thus, these cases could be severely constrained with 
future measurements of the H-polarization mode. 

In such way, we expect that a future statistical analysis 
using latest Planck data, will allow us to constrain the 
values of B for the Newtonian and Wigner’s scheme in 
the case where the collapse occurs when the proper wave¬ 
length of the mode is bigger than the Hubble radius. 

Figures EH 23 and EH show the prediction for the 
temperature fluctuation spectrum for the independent, 
Newtonian and Wigner schemes in the case where k 
TL(r]^). In the independent scheme (Fig. [TUI) , we obtain 
no changes in the CMB temperature power spectrum for 
the different values of B considered and only a small dif¬ 
ference with respect to the fiducial model (in agreement 
with Fig 2]). Therefore, we expect that any value of B 
will explain recent observational data for the values of A 
analyzed in this paper. In the Newtonian and Wigner 
schemes (Figs, fill and fl2l respectively!. for increasing 
value of B we can observe an increase in the value of 
the secondary peaks and a decrease of the value at the 
valleys; the magnitude of the changes depending on the 
value of B, the change is greater for the Wigner scheme. 
Therefore, we expect that a statistical analysis compar¬ 
ing the model predictions with observational data can 
give good constraints on the value of B and in conse¬ 
quence on the time of collapse. In particular, from this 
preliminary analysis, we can expect that the value of B 
should be smaller than 1 Mpc in order to fit recent ob¬ 
servational data. To get a more stringent bound we need 
to perform a statistical analysis using Planck and other 
CMB data and this is left for a forthcoming work [5l j. 
In summary, from the predictions for the CMB tempera¬ 
ture power spectrum shown in this section, we can expect 
that the comparison with recent observational data will 
constrain the values of B for the Newtonian and Wigner 
schemes in both cases analyzed in this paper (i.e. when 
the proper wavelength of the mode is smaller/bigger than 
the Hubble radius at the time of collapse). In contrast, 
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Figure 7. The temperature auto-correlation (TT) power spec¬ 
trum for the independent scheme in the case where k -C 
' H(rIk)- All models are normalized to the maximum of the 
first peak of the fiducial model. The power spectrum of the 
fiducial model is also shown in red. Different values of the col¬ 
lapse time 17 %. = A/k + B are considered, the scalar spectral 
index n s = 0.96; Top: A = — 1CP 1 , Bottom: A = —10 -2 . 


we can expect that any value of B will explain recent ob¬ 
servational data for the independent scheme in the case 
where k H(ri^) and k <C ’H-iVk)- However, we remind 
the reader that this analysis is valid only for the values of 
A considered in this section; a complete analysis studying 
all allowed values of A is in progress [HTJ . 


VIII. SUMMARY AND CONCLUSIONS 

In this paper, we have calculated the primordial power 
spectra for the simplest inflationary model, i.e. a single¬ 
scalar field in the slow-roll approximation, but taking 
into account a collapse of the inflaton wave function (for 
each mode); the motivation for considering an objective 
collapse is to provide a precise explanation for the emer¬ 
gence of an inhomogeneous and anisotropic universe. Un¬ 
like previous works, we have considered a quasi-de Sitter 




Figure 8. The temperature auto-correlation (TT) power spec¬ 
trum for the Newtonian scheme in the case where k -C 'H(rj £). 
All models are normalized to the maximum of the first peak 
of the fiducial model. The power spectrum of the fiducial 
model is also shown in red. Different values of the collapse 
time = A/k + B are considered, the scalar spectral index 
n s = 0.96; Top: A = —10 -1 , Bottom: A = —1CU 3 . 


background to perform our calculations. Even though, 
we have not assumed a particular mechanism for the col¬ 
lapse to happen, we have chosen three different collapse 
schemes, in which the difference between them is due to 
the variable that is affected by the collapse and their cor¬ 
relations. 

The three collapse schemes induced a modification to 
the standard scalar power spectrum of the form P(k) = 
(C/n 2 )k ns ~ 1 Q(\zk\), with n s — 1 = —2 cr + 2<5y and 
Q{\zk\) a function of the time of collapse rj £ (recall the 
definition Zk = kr]f.). 

Moreover, we have characterized the time of collapse as 
r]^. = A/k + B. If B = 0, in all schemes, we can recover 
an equivalent power spectrum that is, for all practical 
purposes, indistinguishable from the standard one. The 
reason for this, is that in this case, Zk = A and the depen¬ 
dence of primordial power spectrum with k is the same 
than in standard inflationary models. The only differ- 
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Figure 9. The temperature auto-correlation (TT) power spec¬ 
trum for the Wigner scheme in the case where k -C 
All models are normalized to the maximum of the first peak 
of the fiducial model. The power spectrum of the fiducial 
model is also shown in red. Different values of the collapse 
time r]% = A/k + B are considered, the scalar spectral index 
n s = 0.96; Top: A = —10 _1 , Bottom: A = —10~ 2 . 


ence in this case is the relation between n s and the slow- 
roll parameters of the inflationary model [see Eqs. (EH 
and m\■ However, we would like to stress that present 
constraints on n s (obtained by comparison with observa¬ 
tional data) that affect the slow-roll parameters can be 
fulfilled by our model as well as the standard one, since 
in both cases the requirement for an efficient inflationary 
stage is ey — Sy 1. 

On the other hand, by assuming B ^ 0, we have found 
small departures from the standard prediction that are 
uniquely determined by the collapse proposal; this is in 
contrast with previous works, in which the introduction 
of the B parameter was used primarily to depart from 
an exact scale invariant spectrum. Moreover, given that 
those works were based on assumptions that led to a 
power spectrum with an spectral index n s = 1, one could 
not tell the precise difference between the dependence on 
k given purely by the dynamics in the aforementioned 




Figure 10. The temperature auto-correlation (TT) power 
spectrum for the independent scheme in the case where k 

All models are normalized to the maximum of the 
first peak of the fiducial model. The power spectrum of the 
fiducial model is also shown in red. Different values of the col¬ 
lapse time 7 ]% = A/k + B are considered, the scalar spectral 
index n 3 = 0.96; Top: A = — 10 2 , Bottom: A = — 10 6 . 

background or by the collapse proposal. 

The primordial power spectrum for the collapse mod¬ 
els obtained in this paper considering B ^ 0 has an addi¬ 
tional dependence on k , which is similar to the one that is 
obtained in standard inflationary models with a running 
spectral index@ However, as it follows from Eqs. (l43l) 
and 65P the dependence on k of the collapse primor¬ 
dial spectrum is different from the standard model one. 
Additionally, as we have mentioned before, the lowest 
multipoles of the temperature_anisotropy are best fitted 
by models with n run ^ 0 jH, |49| . 

We have also shown some plots of the primordial spec¬ 
trum resulting from our schemes for some specific values 


4 The primordial power spectrum with a running spectral index 
has the form: P(k) = A s ( log(k/k 0 ). 
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Figure 11. The temperature auto-correlation (TT) power 
spectrum for the Newtonian scheme in the case where k 
Wir/h). All models are normalized to the maximum of the 
first peak of the fiducial model. The power spectrum of the 
fiducial model is also shown in red. Different values of the col¬ 
lapse time = A/k + B are considered, the scalar spectral 
index n 3 = 0.96; Top: A = —10 2 , Bottom: A = —10 s . 


of the collapse parameters and compared them with the 
standard inflationary model spectrum. We have consid¬ 
ered the case where the associated proper wavelength of 
the modes is bigger and smaller than the Hubble radius 
at the time of collapse. Finally, we have shown the effects 
of introducing a self-induced collapse of the inflaton wave 
function on the CMB temperature fluctuation spectrum. 
For this preliminary analysis, most of the collapse mod¬ 
els proposed in this paper seem to be good candidates to 
explain present data of the CMB fluctuation spectrum. 
In particular, in some cases there are no differences be¬ 
tween the prediction of the collapse models with respect 
to standard inflationary model. However, in other cases, 
for an increasing value of B (affecting the time of col¬ 
lapse) there are important departures from the standard 
model prediction. Therefore, by performing a statisti¬ 
cal analysis using all present observational data from the 
CMB (which is left for future work [5l|) we will be able 




Figure 12. The temperature auto-correlation (TT) power 
spectrum for the Wigner scheme in the case where k 
T-L{r]k)- All models are normalized to the maximum of the 
first peak of the fiducial model. The power spectrum of the 
fiducial model is also shown in red. Different values of the col¬ 
lapse time rg. = A/k + B are considered, the scalar spectral 
index n s = 0.96; Top: A = —10 2 , Bottom: A = —10®. 


to constrain the value of B and thus the acceptable val¬ 
ues for the time of collapse. Finally, we have also shown 
that some particular cases could be discarded without 
performing any statistical analysis. 

We conclude that a more detailed analysis involving 
recent observational data can be used to discriminate 
between the three collapse schemes presented. Never¬ 
theless, from the theoretical point of view, we think that 
the Wigner collapse scheme should be preferred over the 
other two schemes since it takes into account the natu¬ 
ral correlations between the canonical variables that are 
present in the pre-collapse vacuum state. The choice of 
the Wigner distribution to describe these correlations in 
the present setting is justified by some of its standard 
properties regarding the “classical limit,” and, by the 
fact that there is a precise sense in which it is known to 
encode the correlations in question. 
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Appendix A: Deduction of Eq. f[4~f) 

Einstein equations at first order in the perturbations, 
<5Gg = 8 t tGST°, SG° z = 87 tG6T° and SG) = 8nG6Tj, are 
given respectively by 


function gk{y) can be expressed as gkijf) = Re[gfc(? 7 )] + 
ihn{gk{rjj\ (recall that gk = y' k — T~Lyk)- Thus, the real 
and imaginary parts of yk and gk are: 

Re[j/ fe (? 7 )] = (0 7 J v {—krj), (Bl) 

lm[y k {ri)\ = (0 ' ^p^Y v (-ki 7 ), (B2) 

+ y/-ki)J v+ 1 (—fc? 7 ) J, (B3) 


V 2 T - + V) = 47 tG[-<P' 0 2 <S> + $,<$(jf + d^VaH^], 

(Al) 

+ T') = 47rG0i(#,<ty), (A2) 

['£" + H( 24' + $)' + (2 W + U 2 )<f> + iV 2 ($ - 4')]<5j 
- - V) = 4:irG[<j)' 0 8<f)' - - d 0 Ua 2 <^]<5}. 

(A3) 

The perturbations considered 4', < I> correspond to the 
gauge-invariant quantities known as the Bardeen poten¬ 
tials; the perturbation 8(f> and <5<// are the gauge invari¬ 
ant perturbations associated to the inflaton field. For 
the case * / j in Eq. (IA3I) . together with appropri¬ 
ate boundary conditions (more easily seen in the Fourier 
transformed version), leads to 4 1 = 4>; from now on we 
will use this result. 

Equations (lAii) . dm m above, together with Fried¬ 
mann equations, can be manipulated to yield the follow¬ 
ing expression: 


V 2 4t + /i4> = 47tG(o;5(/ + <fi' 0 5(j>'), (A4) 

where n = ’H 2 — H' and oj = 3"H</)g + a 2 d c / > V. Upon use 
of the motion equation for </>□ in the slow-roll approxima¬ 
tion, i.e. 3"Hi/g + a 2 d<f,V ~ 0, implies that u ~ 0. Thus 
m becomes: 


V 2 tf + /i* = 4nG<t>' 0 S^'. (A5) 


+ y/-kr)Y u+1 (-kr})\ (B4) 

where a = 1/2 + v + 1/(1 — e H ). 

With the previous expressions at hand, we can now 

compute ^At)^' 7 (? 7 ^)^ and (r)^ within the in¬ 

dependent and Newtonian collapse schemes. The exact 
form of the quantum uncertainties can be obtained from 
Eqs. (THU) , this is 

(Ayk’ 7 ( 7 7fc))p = ~^\yk{Vk)\ 2 

= ~T,[ j2 M) + Y 2 (\z k \)\ ,(B5) 



L 3 


\9k(rfk) I 5 


L 3 nk 

16 


x 


/ -aJ v (\z k \) 
\ V\Zk\ 



2 


( ~aY v {\zk\) 

V V\Zk\ 



(B 6 ) 


with Zk = hr/fr and 77 ^ the time of collapse for each mode. 
Note that this quantities have the information that the 
background space-time is quasi-de Sitter 

In the case of the Wigner scheme, the explicit formulas 
characterizing the parameters A k(zk) and Ok{zk) are: 


Appendix B: Explicit equations of Sec. inn 
1. Equations of Sec. mm 


A k = (2U) 3/2 '\/® [J 2 (\zk\) + Y 2 (\z k \)] 1/2 [S(N) 


- \ S 2 (\z k \) - 


7 r\Zk\ 


(J3(\ z k\) + Y?(\z k \)) 2 


- 1 - 1/2 


The mode function yk(v) of Eq. (fl2l) . can be expressed 
as Vkiv) = RebfcM] + ilm[j/fc(?;)]; similarly, the mode 


(B7) 
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tan20 fc = [Jl{\z k \) + Y*(\z k \)\ 


S(\z k \)-^(J 2 A\z k \) + Y„\\z k \)y 

x [-2v(Jt(\z k \) + Y?(\z k ))+\z k \ 

X {M\zh\)Jv+i(\z k \) + Y u {\z k \)Y u+1 {\zk\))\, 


where 


(B8) 


F(kr], z k ) = 


kn f ( -aYJ\zk 


M + y/W\Y v+1 (\z k \ 


4 l V VW 

+ y ^^ Jv+1 ^ k ' n ^ 
p/=f- + v/Kp«>(W)) 


(B12) 


S(N)= 1 + ^ \z k \ 2 (Ji{\zk\)+Y v \\z k \)Y 


J„(\zk\) +Y?{\z k \) - \z k \{J v (\z k \)J v+1 (\z k \) 

2> 


Y v {\z k \)Y v+1 {\z k \)) 


(B9) 


2. Equations of Sec. inrci 


In order to deduce Eq. m, we introduce the quantity 
d k ,J = (0|a k -/ |0) that determines the expectation value 
of the field and momentum operator for the mode k at 
all times after the collapse. That is, from Eq. m , we 
have 


(^ J (v))e = VZR-lgk^d*' 1 }, 


(BIO) 


which corresponds to expectation values at any time af¬ 
ter the collapse in the post-collapse state |0). One can 
then relate the value of c? k ’ with the value of the ex¬ 
pectation value of the fields operators at the time of 
collapse {il£ J {vp)e = V^n[y k {r]l)d^ 1 ], <7r k J (v k ))e = 
\f2'R'[9k{Vk)dk’ 1 }- Using the latter relations to express 
d k ’ in terms of the expectation values at the time of 
collapse and substituting it in (IB 101) . we obtain an ex¬ 
pression for the expectation value of the momentum field 
operator in terms of the expectation value at the time of 
collapse 


G(krj, z k ) = T <! J v (\z k \ 




Ujy (I Zk | ) 


- aJ jy ,n + VMMM) 


V-B13) 


Furthermore, given Eqs. (1B11I) and (fTbl) . we can 
find the curvature perturbation in the longitudinal gauge 
within the three collapse schemes. 

For the independent scheme the curvature perturba¬ 
tion is: 


= 


\/L :i TT6vH 


2 5 / 2 M P k 2 

F(k V ,z k )X K1 ^j H [Jl(\z k \) +Y?(\z k \)\ 
G(kij,z k )X K2 Vk 


1/2 


-aJJ\z k 


==~^ + VW\J» + i(\z k \) 


V \/N 

+/ran +1 (i^i)' 


1 1/2 


■ (B14) 


Meanwhile, for the Newtonian scheme we have 




z2 y=f 1 + dMMM) 




+ 'J\zk\Y v+ i(\z k \) 


1/2 

(B15) 


with 


(^k’^))© = F ( k V,z k ){Vk J (.Vk))@ 

+ G{krj,z k ){n FI (r]^))e, (Bll) 


and finally for the Wigner scheme 
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K is (v) = 


e v H 


2 M P k 2 
+ kG(ki k z k ) sin0 fc ]X k , 

with A k = + ix k . 


Afc[F(fc? 7 , Zk) cos 0^ 


(B16) 




ji{kR D )ji{k' Rd) 


k.k' 


fc 3 / 2 fc ' 3 / 2 


x Y^{k)Y lm (k')T(k)T{k') 
x VF(| 2fe |)VF(|^|)^^(^') 3/2 " 


(C3) 


3. Definitions of the functions M(\z k \), 2V(|^fc|) and 
W(\z k \) 


M(\z k \) = ~V\zk\ 
x [J 2 M) 


OiJy 


\ Z k\ 


\/N 

y 2 (\z k \ 


d“ VW\J»+l(\Zk\) 

I 1/2 


(B17) 


with 


C = ( 2 " 11/2r (^ - m\ri? /2 - v ) 2 , (C4) 

where we used the fact that in slow-roll inflation ey — £h- 


Appendix D: On the meaning of the power 
spectrum within the collapse proposal 


IVd^l) = ^/\zk\J v {\zk\) 


X 


/ (| z k 1) 

V V\zk\ 



+ 



1/2 

, (B18) 


In order to obtain the power spectrum in the tradi¬ 
tional inflationary scenario one needs to compute the 
quantum two-point correlation function. That is, if T 
represents the quantum operator associated to the met¬ 
ric perturbation, then the (scalar) power spectrum is ob¬ 
tained from 


W{\z k \) 


27c 1 / 2 

7T 1 /2L3/2 


X 


/ -aJ v {\z k \) 
V \/\zk\ 



+ 


J z>k | Jv (| Zk |)A k sin 0& 


(B19) 


- - 2t r 2 

<0|M k ' |0) = -^-P(k)6(k - k') (Dl) 

On the other hand, let us recall that in general, the 
definition of the power spectrum is given in terms of *I' k , 
i.e. a classical stochastic field not a quantum operator. 
Therefore, the standard approach is based on the identi¬ 
fication 


Appendix C: Equations of Sec. m 


<0|Mv|0)=M k ' (D2) 


The explicit expressions for |a; m |^j L can be found by 
substituting TZ k , given in Eq. (l23l) into Eq. (l28l) . and 
then making the identification |az m |^j L = |aj m | 2 


\air 


12 ind _ I £- 77-2 C ^ ji{kRD)ji{k'R d) 

IML k 3/2 k > 3/2 


x y^(fc)l) m (fc')T(fc)T(fc') (MdzkDMMX^X*, A 
+ N(\z k \)N(\z k ,\)X^X*r?j (fcfc') 3/2 -", (Cl) 


i 12 newt _ i/> 2 C \ ^ Jl{kR D )ji{k Rp) 

I ®lm |ML 1971 2_^ k 3/2 k i3/2 

k.k' 


x^(fc)y Jm (fc')T(fc)T(fc')IV(|^|)A(|^|) 
x X K2 Xl, 2 {kk'f/ 2 ~ v , (C2) 


with 'Ek'I'k' denoting an average over an ensemble of clas¬ 
sical stochastic fields. The justification for the relation 
above relies on arguments based on decoherence and the 
squeezing nature of the evolved vacuum state [H, 0 
(although we do not subscribe to such arguments for the 
reasons exposed in [2, 8|). Nevertheless, the explanation 
about the fact that a single outcome or realization T k 
has emerged is incomplete. 

Therefore, in order to show explicitly the manner in 
which we can obtain an equivalent power spectrum within 
the collapse proposal, we start by focusing on the tem¬ 
perature anisotropies of the CMB observed today on the 
celestial two-sphere and its relation to the scalar metric 
perturbation T, which can be approximated by 

(D3) 

Jo 3 

On the other hand, the observational data are described 
in terms of the coefficients a; m of the multipolar series 
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expansion 


V) = y2 ai m Yim(0, ( f), 

T ° 

aim = J ^(e,ip)Y? m {0,ip)dn, 


(D4) 


here 9 and ip are the coordinates on the celestial two- 
sphere, with Yi m (6,ip) as the spherical harmonics. 

The value for the quantities ai m are then given by 


Aqri l r rj ^ Z- 

aim = —J (2^ MkR D )Y^(k)T(k)^ it (D5) 

with ji(kRjj) as the spherical Bessel function of order 
l and Rd is the comoving radius of the last scattering 
surface. We have explicitly included the modifications 
associated with latetime physics encoded in the trans¬ 
fer functions T(k). The metric perturbation ’I'k is the 
primordial curvature perturbation. 

Now, the problem is that, if we compute the ex¬ 
pectation value of the right-hand side (i.e., identifying 
(4') = ’I') in the vacuum state |0), we obtain 0, while 
it is clear that for any given Z,m, the measured value 
of this quantity is not 0. That is, if we rely in this case 
on the one-point function and the standard identification 
of quantum averages with classical ensemble average, we 
find a large conflict between expectation and observa¬ 
tion. Nevertheless, in the standard approach somehow 
(e.g. by invoking decoherence, squeezing of the vacuum, 
many-world interpretation of quantum mechanics, etc.) 
occurs the transition 'I'k —> ’I'k = ^4e*“ k with a k a ran¬ 
dom phase and A is identified with the quantum uncer¬ 
tainty of i.e. A 2 = (0|'I'k|0), but the random nature 
of ’I'k remains unclear. 

In our approach, the random nature comes directly 
from the stochastic aspects of the quantum dynamical- 
reduction, i.e. from the self-induced collapse. Thus, us¬ 
ing Eq. m in Eq. (ID5I) we obtain 


aim — 


4ni l 

4:TTl l 

3~ 


d 3 k 

en_H_ 
2 M P 


ji(kR D )Y^(k)A(k)^ k 
d 3 k 


= —\l^— I ^ 3 ji(kR D )Y^{k)T{k) 


(ftk) 
k 2 ' 
(D 6 ) 


The previous expression shows how the expectation value 
of the momentum field in the post-collapse state acts as 
a source for the coefficients cq m . 

Furthermore, the angular power spectrum is defined 
by 


Cl = 2[TiEl a ^l 2 - ( D7 ) 

m 

For the reasons presented in Sec. hyi we can identify 
the observed value |a ; m | 2 with the most likely value of 
\ a im\ 2 ML and in turn, assume that the most likely value 
coincides approximately with the average \ai m \ 2 - Thus, 
in our approach, the observed Ci coincides with 

( D8 ) 

m 

From Eq. (ED) we obtain 


«Mf )7 


d 3 kd 3 k' 

(2tt) 6 J A kRD > J A k Rd > 


x Y: m {k)Yi m {k')T{k)T{k')^ k *k' 

2 '■ d 3 kd 3 k' 

— - -Tjr-jl(kRD)jl(k'R D ) 


47T 

~3~ 


(2tt)6 


x Y^YUk'mmk 1 ) 


H (7r k )(TT k >) 


2 M 2 k 2 k' 2 


(D9) 


Consequently using the generic definition of the power 
spectrum [i.e. not relying on the identification (ID2I) ] 


_ 2t r 2 

Mk' = -p-P(fc)«(k - k'), (DIO) 

and also using Eq. (ID91) . the power spectrum, associated 
to 4'k, in our approach is given by 


P ( fc ) = I (v))- (Dll) 

The quantity (Tf k (ri))(n k >( t])) is obtained by using Eq. 
m in the limit —kij —> 0 , i.e. when the proper wave¬ 
length of the modes of interest are bigger than the Hubble 
radius, but taking into account that Eq. (EHD is different 
for each collapse scheme. 
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